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Abstract. 

We present an elementary method to obtain the equations of the shallow-water 
solitary waves in different orders of approximation. The first two of these equations are 
solved to get the shapes and propagation velocities of the corresponding solitary waves. 
The first-order equation is shown to be equivalent to the Korteweg— de Vries (KdV) 
equation, while the second-order equation is solved numerically. The propagation 
velocity found for the solitary waves of the second-order equation coincides with a 
known expression, but it is obtained in a simpler way. By measuring the propagation 
velocity of solitary waves in the laboratory, we demonstrate that the second-order 
theory gives a considerably improved fit to experimental results. 
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1. Introduction 

Solitary waves propagating on the surface of fluids have been extensively studied since 
John Scott Russell discovered them in 1834 p]. They occur naturally in the form of tidal 
bores, but they can be generated in a laboratory as well. The most striking property 
of these solitary waves is their well-distinguished shape which they maintain during 
propagation; they can travel extremely long distances without considerable dispersion 
or dissipation. 

The first theory to successfully explain the existence of solitary waves was developed 
by Korteweg and de Vries [2]. It is based on the shallow- water theory of ideal fluids, 
which predicts linear waves in the small amplitude limit; these waves are slightly 
dispersive. However, nonlinearity is also present due to the convective term appearing in 
the Euler equation. According to the Korteweg— de Vries (KdV) theory, solitary waves 
emerge as a balance between nonlinearity and dispersion. The shape of the waves and 
their propagation velocity can be obtained from the exactly solvable KdV equation. 

The KdV theory can be considered as a first-order approximation to solitary waves. 
Although it explains their unusual properties and agrees well with experiments for 
small amplitude solitary waves, further refinements to the theory are possible. By 
using the systematic expansion method developed by Friedrichs [3], the second-order 
approximation to the solitary waves was found by Laitone |4j. Later on, Fenton extended 
the method to nine orders [5] , while Schwartz reached the 70th-order approximation with 
the aid of computers [6]. 

The derivations of the KdV equation or any higher order approximations are 
absent from many standard textbooks [71 |8], while others giving more complete 
account on the topic use involved mathematical techniques [HE]. In this paper, 
we present an alternative method for treating shallow-water solitary waves. This 
method is mathematically simpler and physically more intuitive; it can be presented 
in any undergraduate course. Based merely on the conditions of incompressibility 
and irrotational flow, we use Bernoulli's law to obtain equations that describe the 
shape of solitary waves in the different order approximations mentioned above. The 
first two equations are solved to find approximate shapes and two expressions for the 
propagation velocities. We test the validity of these expressions by comparing them 
to large amplitude solitary waves in the laboratory, and find that the second-order 
approximation gives much better correlation with the experimental results. 

2. Basic equations 

Let us consider a solitary wave propagating to the left along the —x direction with 
velocity c in an unbounded fluid of ambient height H. Examining it from a co- 
moving system, the stationary flow described in Figure [1] is observed. The origin of 
the coordinate system is placed to the bottom of the fluid under the peak of the solitary 
wave. The horizontal [vx) and vertical {vy) components of the flow velocity are functions 
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Figure 1. An illustration describing the stationary flow of a solitary wave in a co- 
nioving frame. 



of the coordinates x and y, while h{x) denotes the excess fluid height compared to H. 
The amphtude of the sohtary wave is deflned as the value A = h{0). In the limiting 
case of a; — > ±00, it is obvious that Vx{x, y) — > c, Vy{x, y) — > and h{x) 0. 

For many liquids, in particular water, the effect of viscosity can be neglected, 
therefore the problem can be described by the Euler equation. Instead of trying to 
solve it directly, we first recite some important properties of the fiow to be examined: 
first, water is practically incompressible, hence 

Vi; = 0. (1) 

As a consequence, the density p of the liquid is constant. It can be assumed by most 
wave phenomena that the fiow is irrotational so that 

V X i; = 0. (2) 

According to the standard boundary conditions, the vertical velocity must vanish at the 
bottom of the liquid: 

Vy = Q for y = 0. (3) 

It is also clear that the material fiux through the full depth must be the same for all 
vertical cross-sections taken at any x. The constant value of this fiux can be obtained 
by calculating it for the limiting case of x ±00 such that 

,-H+h(x) 

Vx{x,y)dy = cH. (4) 

The last important property of the fiow is that the liquid pressure equals the pressure 
of air on the surface: 

p = Po for y = H + h{x). (5) 

Finding a fiow that satisfies conditions ([I]) — (IS]) is much easier than solving the Euler 
equation directly. 
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3. An iteration scheme 



In this section, we present a method for obtaining an ordinary differential equation for 
h{x) in an iterative sequence of steps. Let us first examine the flow with the velocity 
components 

cH 



and Vy{x, y) = 0. 



(6) 



H + h{x 

It is clear that the conditions ([2]) — (jl]) are fulfllled; however, the divergence of the velocity 
does not vanish. To satisfy condition ([T]), we add a new term to the vertical velocity 
Vy. The new term Avy can be obtained by elementary methods: differentiation of Vx 
with respect to the coordinate x and integration with respect to the single variable y. 
An arbitrary function of x appears after the integration which can be chosen to fulfll 
condition ([3]). The new term obviously does not contribute to the flux, therefore the 
flow with these modifled components satisfles conditions ([1]), ([3]) and (jlj). 

For the new flow, however, the rotation of the velocity does not disappear; to satisfy 
condition ([2]), we now add a new term Avx to the horizontal velocity Vx- The method 
for obtaining Avx is nearly the same as above; after the differentiation of Vy with respect 
to X and the integration with respect to y an arbitrary function of x appears. This can 
be chosen to make the flux resulting from the new term vanish. The vertical velocity 
at the bottom of the water obviously remains zero, hence the flow with the new term 
satisfles the conditions (E]) — (jlj). 







Av(0,Av,) 
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Vx(Av) = 6^0 
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Figure 2. A flow diagram summarizing one iteration cycle, as described in the text. 

At this point, the situation is formally the same as at the beginning: the divergence 
of the flow is non-zero again. It can be assumed, however, that after these steps we are 
closer to the exact flow solution. Once again, a new term must be added to the vertical 
velocity and another one to the horizontal velocity as mentioned above. By following 
this iteration scheme and repeating the steps shown in Figure [21 we obtain inflnitely 
long expressions for both components of the velocity. Assuming that this method is 
convergent, the new terms become less and less signiflcant. The remaining divergence 
and rotation after all the steps therefore approach zero and the limiting flow satisfles 
conditions ©-(HD. 

It is, however, still left to examine whether condition ([5]) is fulfllled. Since viscosity 
can be neglected and the flow is irrotational, Bernoulli's law holds between two points 
of the water surface; one of these points can be chosen to be inflnitely far away, while 
the other is arbitrary with its horizontal coordinate x: 
1 o 1 



Po + pgH + -pc 



Po + pg{H + h{x)) + -p 



H + h{x)) + H + h{x)) . (7) 
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After rearranging this equation, both the air pressure po and the density p of water 
cancel out, hence we obtain 

. 2gh{x) vlix,H + h{x)) , vl{x,H + h{x)) 

^ ~ c^ + c^ • 

This equation determines the shape h{x) of the surface. The velocity components are 
given by the infinitely long expressions obtained from the iteration scheme of Figure [2j 
Equation ([HD is thus an infinitely long ordinary differential equation, which is impossible 
to treat without approximations. 

4. Order of magnitude estimates 

By substituting in principle the infinitely long expressions of Vx{x,H + h{x)) and 
Vy{x, H + h{x)) into equation ([8]) we obtain an equation containing an infinite number of 
complicated terms. These terms can, however, be expanded into series with respect to 
dimensionless quantities much less than unity. Assuming a long solitary wave of small 
amplitude, these quantities are 

and Hh"{x)<^^l, (9) 

for example. The derivatives of h{x) with respect to x are denoted by primes. After the 
expansion the equation is still infinitely long, but its terms are simpler and can more 
easily be classified. Equation ([8]) can formally be written as 

MM_HM,£r„.o, (10) 

n=l 

where all terms denoted by T„ take the similar form: they are the products of the 
function h{x) and its different derivatives multiplied or divided by the appropriate 
power of H to keep them dimensionless. Some of the terms are nonlinear, while those 
containing derivatives are responsible for dispersion. Examples of such terms are 

and H-h"{x). (11) 
Both nonlinearity and dispersion can occur in more complex terms, such as 

■ h\x) ■ h"\x). (12) 

To estimate the magnitudes of these different terms, let us write the shape of the solitary 

wave as 

h{x) = A ■ f{kx) where /(O) = 1. (13) 

The new quantity k is an 'effective wave number', which is inversely proportional to 
the horizontal extension of the solitary wave, while / denotes an appropriately smooth 
unknown function of the dimensionless product kx. The magnitude of the derivatives 
of h{x) can be estimated via equation f|T3l) . For example, the magnitudes of the terms 
in ( ITT]) are 

and H.h\x)^H.Ak\ (14) 
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respectively. It is assumed that the amphtude A is much smaller, while the length ~ l//c 
is much larger than the initial water height H: 

— < 1 and kH <^ 1, (15) 
H 

which are equivalent to the relations in ([9]). 

Solitary waves emerge as a balance between nonlinearity and dispersion; in the 
most simplistic case, it is enough to keep the largest nonlinear and the largest dispersive 
term. These are the terms in ( ITTi) . therefore in this case they must have the same order 
of magnitude: / H"^ ~ HAk"^, i.e. 

4 - k'H\ (16) 

This relation determines the relative magnitudes of the small quantities in fllSp . In the 
following section, we derive the simplest approximation, which corresponds to a case 
where all terms T„ in equation fllUp are neglected except for the terms in (ITT]) . The 
solution can be checked to be in agreement with relation f|T6|) . It is natural to expect 
that the same relation also holds for the more accurate approximations since keeping 
smaller terms does not essentially change the order of magnitude relationships. 



5. First- and second-order equations 

In order to derive the approximate equations, let us return to equation and implement 
the method described in the previous two sections. First of all, we divide the velocity 
components by c, and expand the component Vx into a Taylor series. If we keep terms 
up to the order of A^/H^, we obtain 

, - 1 H ^ and ^ - 0. (17) 

Now we add a new term Avy to the vertical velocity and hence satisfy condition ([1]): 

^ + ^ = 0. (18) 

OX oy 

By keeping terms up to the order of A^ /H^, we obtain 

The arbitrary function g{x) should be chosen to be zero in order to fulfil condition ([3]). 
Next, we add a new term Av^. to the horizontal velocity to cancel out the rotation of 
the flow component ( |T9l) : 

^-^ = 0. (20) 

OX oy 

After the integration with respect to y there remains an arbitrary function g of x: 

Av, _ 1 f 9{^Mdy = - 2h{x))h"{x) - 2h'ixf]y^ + qix). (21) 



c 



dx " 2m 
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The function q{x) should be chosen to make the flux resulting from the component fl21l) 
equal to zero: 

rH+h{x) 1 ^ 2h{x))h"{x) - 2h'{xf\ y^dy + q{x) {H + h{x)) = 0. (22) 

Hence the velocity components after the first iteration cycle are given by 

h{xY^ 



^=(1 



and 



h{x) ^ h{xY 



H 



2H^ 



{H - 2h{x))h"{x) - 2h\xf 



\{H + h{x)f 



h'ix)f^ 2hix)\ 



(23) 
(24) 



The substitution of these components into equation ([8]) yields 



2gh{x) 



+ ^ 1 + 



h{x) h{xY h{xY 
h{x) 



H 



h'{x) 



(25) 



[H - 2h{x))h"{x) - 2h\xf 

if we keep terms up to the order of A^/H^. The magnitude of the different terms can 
be estimated via equations f|T3l) and f|T6|) . similarly as for the terms in ([1] 
h{x) A 



and 



H 

h{xf 



-Jp- -H-h ix)r^ — 



h{x) ■ h'\x) ~ h\xf 



(26) 



(27) 



Rearranging the right-hand side of equation (l25l) and keeping the terms up to the order 
of A^/H^ leads to the form ([ID]): 

2gh{x) 2h{x) 



where 



and 



H 

3h{xf 2 
3 

4/i(x^^ 



+ Ti + T2 = 0, 



H ■ h"{x) 



'^hix) ■ h"ix) - -h'ixf 



are on the order of A^ /H"^ and A^^/H^, respectively. 

The simplest approximation is to neglect both Ti and T2 in equation 
case, we obtain the linear equation 



2g 



H 



h{x) = 0, 



(28) 

(29) 

(30) 
In this 

(31) 



which describes the shallow-water linear waves of arbitrary shape. These non-dispersive 
waves propagate with the well-known velocity of 

c = Co = J^, (32) 
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which converts equation (1311) into an identity for any h{x). 

Nonhnear and dispersive terms are, however, both necessary for a sohtary wave 
solution. The simplest approximation describing such solutions can be obtained by 
keeping Ti ~ /H"^ and still neglecting T2 ~ A^/H^ in equation fl28|) : 
2gh{x) 2h{x) 



H 



0. 



(33) 



The first-order approximate equation (|33|) is, to the given approximation, equivalent to 
the KdV equation, which is generally used to describe solitary waves |2J. Since relation 
(152]) is a zeroth-order solution for the propagation velocity c, it is worth introducing the 
small dimensionless quantity 
1 3H 



< 1 



(34) 



and approximating the first two terms in equation (1331) as 



(2_gH_ 

V C2 



h{x) 



4-4 1 + 



h{x) 



4- 



h{x) 



4- 



4c \ h{x) 



.(35) 



The approximation (1 + e)° ^ 1 + ae appears in the second step, which is valid to first 
order in e. After differentiation and multiplication by H/A we obtain 

c 



h'ix) - 



j=h\x) + ^h{x) ■ h'ix) + \H^h"'( 



X) 







(36) 



from equation fl33|) . This ordinary differential equation is then transformed into a partial 
one with the independent variables x and t. Derivatives in all terms of equation (I5B1) 
must be substituted with a partial derivative either with respect to x or with respect 
to t. In the latter case, however, the term also must be divided by c so that equation 
(1361) can be restored by seeking the solution in the special form h{x + ct). Since the 
propagation velocity c is different for every special wave solution, it must not explicitly 
appear in an equation describing such a wide range of phenomena. Only the second 
term in equation (l36l) should therefore contain a time derivative. By leaving spatial 
derivatives in all other terms we obtain 
dh 1 

dx 



dh 3 ^ 9/i 1^2^^^ 
dt 2H dx 6 (9x^ 



0, 



(37) 



the KdV equation describing waves propagating to the left. 

It can be verified that further steps after the first iteration cycle do not result in new 
terms with magnitudes up to the order of A^ /H^. For the second-order approximation, 
however, we must collect all terms with magnitudes up to A^/H^; this requires one more 
iteration cycle. If we keep terms up to the order of A^/H^, the new term Avy added to 
the vertical velocity is given by 



Av' 



diAv.,] 



dy 



h"'{x) y^ -y{H + h{x)f + g{x). (38) 



dx " 6H 

The arbitrary function g{x) is once again chosen to be zero so that the flow fulfils 
condition ([3]). By still keeping terms up to A^/H^, the new term Av'^ added to the 
horizontal velocity reads as 



Av' 



d{Av'y[ 
dx 



■dy 



1 

'2AH' 



h^^\x) \y^ - 2y^{H + h{x)f] + r(x), (39) 
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ISy^ - 30|/2(i/ + h{x)f + 7{H + h{x)f . (41) 



where the function r(x) is chosen to make the flux resulting from the flow component 
fl39|l vanish: 

rH+h{x) l , , r on 

/ -^^^'^(^) \y' - + f'(^)) ]dy + r{x){H + h{x)) = 0.(40) 

Jy=0 2AH L -I 

Hence the new term Av'^ is given by 

A< _ /iW(x) 
^ ~ ~ 360iJ 

Similarly as after the first iteration cycle, the velocity components of the net flow 
obtained from expressions fl23l) . fl2^ . fl38|) and PTI) are substituted into equation ([8]), 
which is then rearranged to the form flTU]) . By keeping terms up to the order of A^/H^, 
the three terms of T2 in equation (!30|) appear with respect to equation (!33|) . One new 
term, 2H^h^'^\x)/A5, also results from the component (14T]) of the second iteration cycle. 
The second-order approximate equation thus takes the form 

^+"^+3^'^^"^^ ^-3^(a;)-^ (a;)-3^ (a;) +^^'^^ Ha;) = 0.(42) 

By continuing the same method it is possible to obtain more accurate (higher order) 
approximations without too much difficulty. The calculations, however, become far more 
complicated, hence these further equations are beyond the scope of this paper. 

6. Solitary wave solutions 

In this section, we first recover the well-known properties of the KdV solitary waves 
from equation fl55]) . then find a solitary wave solution for the second-order approximate 
equation (H2l) . By using abbreviation (IMI) . equation (133!) reads as 

h(x) h(x)'^ 2 ,,,, , , , 

- 2e-^ + 3-^ + -H ■ h"{x) = 0. (43) 

After a multiplication by h'{x) it can be integrated to give 

- + + ■ h\xf = B. (44) 

In solitary waves h'{x) for h{x) 0, hence the integration constant B is zero. By 
expressing h'{x) from equation fj44|) we obtain the separable differential equation 



h'{x) = ±J—-h{x)^eH-h{x). (45) 



The solution for equation fl45l) can be obtained by elementary methods: 

h{x) = eH -cosh-^ {k{x -C)), (46) 
where k is the 'effective wave number' given by 
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Figure 3. Shape of the sohtary wave of the first-order equation ([33ll with water height 
H = 1 and amplitude A = 0.1. 



The integration constant C is only responsible for shifting the wave along the x-axis, 
therefore it can be chosen to be zero without any restriction. By comparing equations 
([13]), die]) and (HZD we obtain 



(48) 



between the amplitude and the length of the solitary wave, clearly verifying the validity 
of relation (fT6|) . The propagation velocity c can be expressed from equations (IMIl and 
by means of the following approximation: 



A 



2h)' 



(49) 



'1-A/H 

Relations fHHj) and P9l) are well known from the KdV theory of solitary waves [2], as 
well as the shape ( 146|) of the wave, as plotted in Figure [3l 

Unlike the previous one, the second-order approximate equation (H2l) can not be 
solved in such an elementary way. We first examine the function h{x) asymptotically, 
i.e. in the region where h{x) —>■ 0. By neglecting nonlinear terms in equation P2|) we 
obtain 



2e^ + lH.h"ix) + ^H'h^'\x) 



0, 



(50) 



with abbreviation used once again. Seeking the solution in the exponential form 
exp(fcx) gives four distinct roots: 

1/2 

(51) 



ki2 = ±— 

1,2 ^ 



■ -1 



and 



k 



3,4 



±- 



1 



/15 

'H\ 2 ^ , 
The magnitude of these roots is estimated for e ^ 1 to be 



1/2 



ki2 ~ ±— v3e and 
H 



'-3,4 



±-Vl5. 



(52) 



(53) 



While the first two roots are appropriately small, the last two describe sinusoidal waves 
with wavelength comparable to the ambient water height H. Accepting such waves 
would contradict the second relation of (IT^ , which is an assumption leading to equation 



Higher order corrections for shallow-water solitary waves 11 



A 



h(x) 









Figure 4. Definition of the quantities h{xo) and A, and the numerical solution of 
equation (glD with xq = -6.78, h{xo) = 0.01 and e = 0.1. 



( 142|) . The solitary wave solution therefore must approach the real exponential functions 
corresponding to ki^2 as h{x) — > 0: an increasing one {ki > 0) for x — >• — oo and a 
decreasing one < 0) for x +oo. 

In the region between the two limiting cases, we solve equation f j42l) numerically. 
The initial values for the function h{x) and its derivatives can be set according to the 
asymptotic form described above. If we start from a small value h{xo) located at a large 
negative xq, the derivatives have initial values 

/iW(a;o) = /i(a;o)-fc^ (54) 

If the initial derivatives are not set according to this, short sinusoidal waves appear, 
leading to the above-mentioned contradiction. As long as /i(xo) is sufficiently small, 
changing its value only leads to a shift along the x-axis. For every e therefore we can 
find a solitary wave solution with a unique shape and amplitude. The latter is read off 
as the maximum value of the numerically obtained function h{x), as shown in Figure 
m The results obtained for different e values are summarized in Table [H along with 
two simple functional forms of A/H, the second of which appears to be a really good 
approximation of e. 



e 


A/H 


A/H - A^/H'^ 


A/H -21A^/20H^ 


0.01 


0.0101072 


0.010005 


0.0099999 


0.02 


0.0204382 


0.020020 


0.0199996 


0.03 


0.0310081 


0.030047 


0.0299985 


0.04 


0.0418337 


0.040084 


0.0399961 


0.05 


0.0529342 


0.050132 


0.0499921 



Table 1. The quantity e versus different functional forms of A/H obtained from 
numerical solutions of equation 



The relationship between e and the amplitude A indicates that a correction should 
be added to relation (HHI) in the case of second-order solitary waves: 
A 21A^ 



^ H ^^^^ 
Although Tabled] only shows values of e much smaller than unity, relation (!55|) remains 
approximately valid for larger values. After the critical value of e ~ 0.25, however, there 
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Figure 5. Solitary wave of tlie second-order equation (|42|) obtained numerically with 
H = 1 and A = 0.2 (a); the same wave plotted together with the analytically known 
first-order wave marked by dashed line (b). 



are no solitary wave solutions; the function h{x) diverges exponentially. The highest 
possible solitary waves are found to have the ratio A/H 0.5. While equation (l33ll 
formally allows solitary waves of arbitrary height, equation (j^2|) contains a hint on the 
experimentally observed instability (wave breaking) of high solitary waves. 

The propagation velocity can once again be expressed from equations and 
( l55ll . This time, however, we have to keep terms up to the order of A^/H^ to obtain 
the velocity for the solitary waves of the second-order equation ( H2l) : 

This result is equivalent to that obtained by Laitone with a systematic expansion of 
the flow components [1]. Note that the velocity given by expression fl56|) is smaller 
than the velocity in (H9l) : the maximum correction, belonging to A/H ^ 0.5, is about 
(3/80)/(5/4) PS 3%. The new terms in equation (H2!) with respect to equation (133!) are 
only small corrections, thus the new solitary waves plotted in Figure [5]^a) are similar 
to those of the KdV theory. Comparison in Figure [5](b) indicates that the second-order 
solitary waves are slightly longer. 

To summarize, the results obtained via the numerical methods show that the 
solitary waves of the second-order equation (H2!) are slower and longer than the 
corresponding KdV waves. 

7. Velocity measurements 

In order to experimentally examine the validity of the theory, we measured the 
propagation velocity of large amplitude solitary waves and compared the results with 
relations (l49l) and (!56l) . Solitary waves propagating in a long, narrow glass tank were 
observed with a CCD camera standing in the direction perpendicular to the propagation. 
The experiments were carried out with coloured tap water to make the contrast between 
the fluid and the background larger. The rightmost part of the tank shown in Figure [6] 
was separated by a lock-gate and contained water of height H', larger than the height 
H of the ambient fluid. After pulling out the lock-gate, this difference of water height 
generated a solitary wave propagating to the left. Solitary waves of different amplitude 
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were produced by changing the modified water height H' . The camera was placed to 
look at a region 6 m away from the right end of the tank to avoid transient phenomena. 
Reflected waves propagating to the right were also recorded. 



L 



d 

< > 



V 



H 



H' 



Figure 6. The initial set-up used to produce solitary waves (L = 10 m, d = 13 cm, 
H = 10.3 cm, while H' varied between 17 cm and 24 cm). 



The pictures taken by the camera were afterwards digitally evaluated by a computer 
program, which calculated the displacement of the solitary wave between subsequent 
images. The profile of the wave could easily be obtained by finding the border separating 
the dark and light regions of the picture. Intersections of the surface and different 
horizontal lines were used to monitor the solitary wave on subsequent images, as shown 
in Figure [71 

The exact value of the water height H and the amplitude A could also be read off 
by the computer program. Both values are obtained as the difference of appropriate 
vertical coordinates: the former is the distance between the bottom of the tank and the 
surface far in front of the solitary wave, while the latter is the maximum value of the 
surface with respect to H. In the case of direct solitary waves, both values H and A are 
known with the accuracy of ±1.5 mm. When examining reflected waves, the maximum 
value used to calculate the amplitude A has a larger deviation on subsequent images, 
hence the error of A increases to ±2.5 mm. 



Figure 7. Two photographs taken in rapid succession {At — 0.0635 s) with arrows 
indicating the methodology for determining the velocity of the wave: the displacement 
Aa; = xi — X2 = 7.54 cm gives c = 119 cm s~^. 

Since the time interval between two pictures is known with a negligible relative 
error of ~ 0.1%, the velocity can easily be obtained from the displacement. We used 
several subsequent images to flt a line on the relative displacements to reduce the error 
in the velocity. This error was estimated as the standard deviation of the different 
velocities obtained from tracing the left and right intersections of 9-26 horizontal lines, 
their number depending on the amplitude of the solitary wave. The average values of 
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Figure 8. Dimensionless propagation velocity of the solitary waves as a function of 
A/H; the continuous lines mark the theoretical predictions and ([55]) . 



these velocities for the same wave are plotted non-dimensionally in Figure [HI with the 
vertical error of the points coming from two main sources. Besides the deviation in 
the velocity, the error in the water height H must also be taken into account when 
calculating the Froude number 

Fr = -£=. (57) 

VgH 

The theoretical predictions (1491) and (1561) are also plotted for comparison. Figure [S] 
clearly shows the coincidence of the measured points and relation (l56l) . and indicates 
therefore the validity of the underlying second-order theory. 



8. Conclusions 



The iteration scheme presented in Section [3] can be used to obtain equations describing 
solitary waves in arbitrary order of approximation. These equations are ordinary 
differential equations, therefore they only give solutions of stationary shape (one-soliton 
solutions). By implementing the elementary iteration scheme of Section [3], we obtained 
the first- and second-order approximate equations. 

The first-order equation was solved in Section [6] to recover the well-known solitary 
waves of the KdV theory. By finding the solitary wave solutions of the second- 
order equation and comparing them to those of the KdV theory, we found that the 
second-order solitary waves are slightly slower and longer than the corresponding KdV 
waves. The expression obtained for their propagation velocity is equivalent to the result 
derived by Laitone Velocity measurements of real solitary waves show a really 
good correlation between the experimental results and the second-order theory, which 
is therefore experimentally verified to be valid. 
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It is worth mentioning that further refinements to the theory do not give observable 
corrections to the propagation velocity. By taking the third-order approximation we 
would obtain c/ ^gH = 1 + A/2H - 3A^/20H^ + ?,A^ /hQH^ In the case of the 
highest solitary waves with A/ H ^ 0.5, the last term is only responsible for a maximum 
correction of ~ 0.5%, which is very difficult to demonstrate in an experimentally reliable 
way. The second-order approximation is sufficient to describe shallow-water solitary 
waves in all realistic situations. 

Acknowledgments 

We would like to thank Imre M Janosi, J W A Robinson and Tamas Tel for helpful 
discussions, and Balazs Gyiire for technical assistance. This work was supported by the 
Hungarian Science Foundation (OTKA) under grant NK72037. 

References 

[1] Drazin P G and Johnson R S 1989 Solitons: An Introduction (Cambridge University Press) 

[2] Korteweg D J and de Vries G 1895 Phil. Mag. {Ser. 5) 39 422 

[3] Friedrichs K O 1948 Commun. Pure Appl. Math. 1 81 

[4] Laitone E V 1960 J. Fluid Mech. 9 430 

[5] Fenton J 1972 J. Fluid Mech. 53 257 

[6] Schwartz L W 1974 J. Fluid Mech. 62 553 

[7] Landau L D and Lifshitz E M 1987 Fluid Dynamics (Oxford: Pergamon) 

[8] Kundu P K, Cohen I M and Hu H H 2004 Fluid Mechanics (New York: Academic) 

[9] Ludu A 2007 Nonlinear Waves and Solitons on Contours and Closed Surfaces (Berlin: Springer) 



